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ABSTRACT 


This  report  focuses  on  the  use  of  certain  detection  statistics, 
in  particular  the  instantaneous  amplitude  or  envelope,  and  the  instantaneous 
power,  in  the  design  of  a controlled  false  alarm  rate  detector.  To  achieve 
false  alarm  rate  control,  the  detection  statistic  must  be  stationary,  but  need 
not  be  Gaussian.  Parameters  of  a Gaussian  process  can  be  conveniently  trans- 
formed into  a stationary,  normalized  detection  statistic.  For  an  envelope  de- 
tector the  normalization  consists  of  dividing  the  envelope  by  a long-term  (1  to 
2 minutes)  estimate  of  the  RMS  noise  level;  the  instantaneous  power  is  normal- 
ized by  dividing  it  by  the  long-term  average  power. 

The  stationarity  of  the  normalized  envelope  is  demonstrated 
with  one  hour  of  uninterrupted,  Korean  Seismic  Research  Station  single- site, 
short-period  noise,  by  compiling  distributions  of  the  base-ten  logarithm  of  the 
normalized  envelope.  The  envelopes  closely  follow  the  Rayleigh  distribution 
which  is  the  theoretical  distribution  for  envelopes  of  a Gaussian  process.  Re- 
taining this  model  at  the  5%  significance  level  establishes  adequate  control  of 
the  false  alarm  rate.  Typically,  an  envelope  threshold  of  12  dB  above  the  RMS 
noise  level  should  yield  13  false  alarms  per  hour  for  a single  detection  trial. 
The  use  of  multiple  detection  trials  and  multiple  detection  criteria  should  re- 
sult in  a lower  false  alarm  rate. 

Various  detector  design  optimization  approaches  are  discussed. 
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SECTION  I 
INTRODUCTION 


In  the  concept  of  a world-wide  seismic  surveillance  system  such 
as  described  by  Sax  et  al.  (1974),  reliable  automatic  detection  and  timing  of 
seismic  event  signals  are  essential.  If  the  system  is  to  be  optimized  for  event 
detection  then  the  effects  of  missed  signals  and  false  alarm  detections  must  be 
taken  into  account.  The  false  alarm  rate  (FAR)  must  be  made  controllable, 
i.  e.  , it  must  be  made  constant  and  predictable  for  a given  threshold  setting 
of  the  detection  statistic,  and  should  be  as  low  as  possible  for  a given  proba- 
bility of  detection.  Inversely,  the  probability  of  detection  should  be  as  high 
as  possible  for  a given  FAR.  This  is  a general  requirement  of  any  communi- 
cation system  (Schwartz  et  al.  , 1966;  Van  Trees,  1968).  This  requirement  is 
also  applicable  to  the  seismic  surveillance  problem  in  that  P-wave  transmis- 
sions from  a seismic  event  convey  the  timing  information  needed  to  locate  and 
time  transient  seismic  events,  to  edit  time  windows  containing  signals  assoc- 
iated with  the  event,  and  finally  to  identify  the  source  as  an  earthquake  or  as 
an  explosion.  For  a seismic  surveillance  system,  an  example  of  optimum 
threshold  control  was  worked  out  in  a study  on  feedback  and  parameter  up- 
date aspects  in  such  a system  (Unger  et  al. , 1974). 

Several  recent  developments  were  aimed  at  providing  front-end 
short-period  (SP)  signal  detection  algorithms  which  could  satisfy  the  above 
mentioned  criteria.  Swindell  and  Snell  (1977)  designed  an  automatic  power 
detector  with  a controllable  FAR,  by  transforming  the  detector  output  into  a 
zero-mean,  unit- variance  normal  random  variable  (r.  v.  ).  At  the  Korean 
Seismic  Research  Station  (KSRS)  array,  a comparison  was  made  between 
automatic  signal  detections  and  an  analyst's  detections.  The  results  indicated 


• + 


that  90%  of  the  analyst's  signal  picks  were  detected  if  the  threshold  of  the 
automatic  detector  was  set  at  8.4  dB  signal-to-noise  ratio  (SNR)  and  50%  of 
the  analyst's  picks  were  detected  with  a threshold  set  at  9.  6 dB.  These  fig- 
ures are  based  on  a fixed  variance  of  the  noise  power  of  3 dB.  The  corres- 
ponding FAR  for  these  two  threshold  settingsf  when  employing  eight  beams,  is 
25  false  alarms  per  hour  (FA/H)  at  the  8.4  dB  threshold  setting  and  4 FA/H 
at  the  9.  6 dB  threshold  setting. 

The  receiver  operating  characteristics  (ROC),  i.  e. , the  prob- 
ability of  detection  versus  the  probability  of  false  alarm  (or  the  FAR)  char- 
acteristics, can  be  improved  by  processors  which  enhance  the  SNR  of  the  de- 
tector input  waveform.  An  example  of  such  a processor  for  array  waveforms 
is  a recently  improved  adaptive  beamforming  (ABF)  algorithm  (Shen,  1977, 
1978)  with  typical  SNR  gains  of  5 to  10  dB.  This  algorithm  can  be  implemented 
in  front-end  processors  such  as  the  station  processor.  Other  methods  of  SNR 
enhancement,  of  course,  are  various  forms  of  pre-filtering  such  as  bandpass 
or  Wiener  filtering.  The  problem  with  these  latter  methods  is,  however,  that 
if  designed  to  apply  to  signals  of  diverse  spectral  contents,  their  processing 
gain  cannot  be  very  high. 

Another,  independent  effort  (Unger,  1978)  evaluated  the  use  of 
the  instantaneous  amplitude  or  envelope,  the  instantaneous  phase,  and  the  in- 
stantaneous frequency  in  the  automatic  detection,  timing  and  measurement  of 
seismic  signals.  This  study  resulted  in  the  design  of  an  automatic  phase  de- 
tector and  timer  for  long-period  (LP)  surface  waves,  and  of  an  automatic  en- 
velope detector  and  timer  for  SP  bodywaves.  The  latter  was  designed  to 
closely  follow  an  analyst's  signal  detection  logic.  An  initial  evaluation  of  this 
detector  showed  promising  results:  for  a given,  but  not  necessarily  typical, 
set  of  Norwegian  Seismic  Array  (NORSAR)  single- site  recordings  of  Eurasian 
earthquake  and  presumed  explosion  waveforms,  the  SP  envelope  detector, 
operating  at  a peak- signal- to-peak- noise  ratio  threshold  of  3 dB,  detected  69% 
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j of  16  signals  having  less  than  12  dB  SNR,  with  a FAR  of  7 FA/H;  the  RMS 

timing  error  was  about  0.  2 sec.  Remaining  technical  problem  areas  for  this 
detector  are  detecting  and  timing  emergent  signals  and  signals  interfering 
1 with  the  coda  of  a preceding  signal,  and  controlling  the  FAR.  The  solution  of 

i these  problems  must  precede  any  extensive  evaluation  of  the  detector  to  deter- 

■ mine  realistically  its  full  operating  characteristics  and  detection  capability. 

| This  report  focuses  on  one  of  the  above  problems,  the  FAR 

control  of  the  envelope  detector.  The  study  analytically  and  empirically  de- 
scribes the  probability  density  function  of  envelope  measurements  of  seismic 
noise.  If  the  model  fits  the  data  sufficiently  close,  it  can  be  used  to  establish 
algorithms  for  the  efficient  control  of  the  FAR. 

The  report  is  organized  as  follows.  Section  II  contains  our 
theoretical  development.  We  first  show  briefly  how  FAR  control  is  established 
in  principle  for  a given  detection  statistic  distribution.  We  then  derive  norm- 
alized detection  statistics  and  their  distributions  for  a Gaussian  noise  process,  ] 

and  describe  how  these  enable  FAR  control.  In  Section  III  we  present  empir- 
ically gathered  distributions  and  observations,  and  test  the  goodness-of-fit  of 
our  theoretical  models.  In  Section  IV,  we  interpret  this  information  in  terms 
of  detector  design  considerations.  The  study  is  summarized  in  Section  V and 
related  literature  is  listed  in  Section  VI. 

1 ?: 
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SECTION  II 
THEORY 


A.  INTRODUCTION 

In  this  section,  we  describe  the  control  of  a detector's  false 
alarm  rate  (FAR)  for  Gaussian  noise.  First,  we  describe  the  general  pro- 
blem of  maintaining,  in  an  average  sense,  a constant  false  alarm  rate  (CFAR) 
for  a signal  detector.  Next,  we  discuss  the  problem  of  normalizing  a time- 
varying  Gaussian  process  to  a zero-mean,  unit- variance  or  N(0,  1)  normal 
process.  We  then  derive  the  normalized  distribution  of  two  signal  detection 
parameters  and  their  base-ten  logarithms:  the  instantaneous  amplitude  or 
envelope  and  the  instantaneous  power.  In  each  case  we  show  how  these  norm- 
alized random  variables  (r.  v.  ) may  be  used  in  the  control  of  the  FAR  of  a 
seismic  signal  detector,  for  a given  noise  frequency  band. 

For  an  excellent  summary  of  related  tutorial  material  and  an 
extensive  list  of  literature,  we  refer  to  the  first  chapter  of  Schwartz  et  al. 
(1966).  This  material  is  the  basis  for  our  development. 


B.  FALSE  ALARM  RATE  CONTROL 

If  the  noise  distribution  of  a detection  statistic,  x,  is  known, 
the  FAR  of  a detector  can  be  controlled  by  setting  a threshold  xT<  For  values 
*>XT  the  noise  hypothesis  is  rejected  and  a signal  detection  is  declared.  The 
probability  of  a false  alarm  is 


J p(x)dx. 


(n-i) 


where  p(x)  is  the  probability  density  function  of  the  r.  v.  , x,  in  the  absence 
of  a signal.  The  false  alarm  rate  then  is 


FAR  = 


P^.  3600 
F 


FA/H 


(II- 2) 


where  dT  is  the  interval  fin  seconds)  between  independent  detection  trials. 

This  interval  depends  on  the  noise  frequency  band  (Sections  III  and  IV).  For 

our  study  we  assume  dT  = 0.  36  sec;  in  Equation  (II-2)  the  FAR  then  cor- 
4 

responds  to  10  P . 

F 

If  a constant  FAR  is  required  to  detect  signals  in  non- station- 
ary noise,  the  threshold  x^,  must  adapt  to  changes  in  the  parameters  of  the 
noise  probability  density  function  and  therefore,  the  threshold  will  have  to 
change  with  time.  In  contrast,  a stationary  detection  statistic  yields  a 
constant  FAR  for  a given  threshold  setting.  In  many  cases,  a stationary  de- 
tection statistic  can  be  obtained  by  means  of  a suitable  transformation.  It  is 
also  convenient  to  express  the  threshold  settings  logarithmically,  for  instance 
in  dB  above  the  RMS  noise  value.  These  transformations  are  discussed  in  the 
following  subsections,  for  a Gaussian  r.  v. , and  for  the  envelope  and  the  in- 
stantaneous power  of  a Gaussian  process. 

C.  NORMALIZED  GAUSSIAN  NOISE  PROCESS 

If  noise  stems  from  a large  number  of  independent  (or  non- 
interacting) particles  or  sources,  the  noise  process  tends  to  be  Gaussian,  i.  e.  , 
the  parameters  involved  in  this  process  are  r.  v.  whose  variation  is  ruled  by 
the  Gaussian  probability  density  function: 


P(*)  = 


V~2na 


(II- 3) 


1 


where 


x is  a r.  v. 

4 is  the  mean  or  average  value  of  x 
<T  is  the  standard  deviation  of  x. 

Examples  of  a Gaussian  process  are  thermal  noise,  shot  noise.  Brownian  mot- 
ion. In  seismology,  if  instantaneous  random  pressure  changes  in  the  earth  are 
caused  by  a large  number  of  independent  sources,  the  resulting  seismograms 
show  waveforms  in  which  the  instantaneous  value  n(t)  is  a Gaussian  r.  v.  Re- 
sults of  this  nature  were  found,  among  others,  by  Sax  (1965);  this  article  also 
lists  related  literature.  The  variance  of  this  r.  v.  equals  the  total  power  N, 

or,  equivalently,  the  square  of  the  RMS  level  of  zero-mean  noise  (Schwartz 
et  al.  , 1966): 


RMS2  . 


(II- 4) 


In  Equation  (II- 3)  both  U and  cr  may  vary  with  time.  In  noise  seismograms  the 
mean  is  usually  zero;  the  total  noise  power  in  general  varies  with  time. 

In  designing  a signal  detector  with  a controllable  false  alarm 
rate  it  is  desirable  to  transform  the  detection  statistic  into  a r.v.  which  is 
stationary,  i.  e.  , whose  statistics  do  not  vary  with  time.  In  a Gaussian  pro- 
cess this  can  be  achieved  by  performing  the  transformation  (Lacoss,  1972); 

z - x ~ * 

o ‘ (II-  5) 

This  changes  Equation  (II- 3)  into: 


p(z)  = i- 

y/Fn 


(H-6) 


H-3 


which  always  has  zero  mean  and  unit  variance,  independent  of  time. 

This  z -statistic  was  adopted  by  Sax  et  al.  (1974)  as  a preferred 
detection  statistic  for  automatic  front-end  detection  in  a world- wide  seismic 
surveillance  system.  Swindell  and  Snell  (1977)  implemented  this  concept  to 
control  the  FAR  of  an  automatic  power  detector.  This  was  accomplished  by 
transforming  the  time- varying,  near-Gaussian  distributed  logarithm  of  the 
short-term  average  power  (log  STA)  in  this  manner  into  a stationary  detection 
statistic. 


D.  FAR  CONTROL  WITH  A NORMALIZED  ENVELOPE 


1.  Envelope  Distribution 


It  can  be  shown  (Schwartz  et  al.  , 1966)  that  the  envelope  R(t) 
of  a Gaussian,  zero-mean  random  process  n(t)  has  a Rayleigh  probability 
density  function: 


P(R) 


(H-7) 


where  N = is  the  total  noise  power  over  the  noise  frequency  band.  This 
function  is  sketched  in  Figure  II- 1. 


below: 


Some  properties  of  the  Rayleigh  distribution  are  summarized 


(i)  The  entire  distribution  is  characterized  by  a single  parameter, 
i.  e. , the  most  likely  value  Rj  . This  is  the  value  for  which 
p(R)  is  maximum.  By  differentiation  we  obtain: 


Ri  = V'N  = RMS 


(H- 8) 


II-4 


FIGURE  II- 1 

RAYLEIGH  DISTRIBUTION 
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pmJR)  = p<v  = e 


-i/(2V5T) 


(II-  9) 


(ii)  The  median  or  50%  cumulative  distribution  point  occurs  at 


R , = 1.185VN  . 

med 


(n-10) 


(iii)  The  mean  or  first  moment  is 


li  = V“rN  = 1,25v^ 


(iv)  The  second  moment  is 


m = 2N 
2 


(v)  The  variance  is 


2 2 IT  . 

°R  = m2  ' mi  = N(2  ‘ ‘ 


(II- 11) 


(n-i2) 


(n-13) 


From  Equations  (II— 11)  and  (11-13)  follows 


0. 525  , 


(II-14) 


independent  of  the  noise  power.  This  value  may  be  used  to  test  if  the  envelope 
of  a given  waveform  is  indeed  Rayleigh-distributed. 

2.  Normalized  Envelope  Distribution 

To  enable  convenient  FAR  control  of  an  envelope  detector  we 
proceed  to  normalize  the  Rayleigh  distribution  by  making  it  independent  of 
the  noise  power.  For  this  we  substitute 


R _ R 
y/W  RMS 


(H- 15) 


. r*  . 


iW  *!*•'****••* 


in  Equation  (II- 7),  and  require 


p(R)dR  = p(u)du  , 


(H- 1 6) 


i.  e. , the  probability  that  R has  a certain  value  must  be  the  same  as  the  pro- 
bability that  u has  the  corresponding  value.  This  yields 


p(u)  = ue 


(H-17) 


This  function  is  sketched  in  Figure  II-2.  We  observe  that  the  most  likely 
value  of  u always  equals  one: 


uj  = 1 


independent  of  the  noise  power,  and  that 


(11-18) 


Pmax{u)  = P(V  = 


(II-19) 


3.  Normalized  Log-Envelope  Distribution 

Since  a detector  is  usually  operated  in  a logarithmic  fashion, 
i.  e.  , setting  thresholds  in  terms  of  dB  SNR,  we  now  derive  the  normalized 
log-envelope  distribution.  For  this,  w’e  substitute 


y = log  u = log  R - log  V N = log(R/RMS) 


(11-20) 


so  that 


d_y  _ k 
du  u 


(11-21) 
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ur 

FIGURE  II- 2 

NORMALIZED  RAYLEIGH  DISTRIBUTION 
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where 


k = log  e = 0.43429  . 


(11-22) 


Requiring 


we  obtain 


p(y  )dy  = p(u)du 


"k"p(u) 


(H-23) 


(H-24) 


where 


01-25) 


- 1 y 

u s antilog(y)  — log  y = 10  . 


(11-26) 


The  maxima  and  minima  of  p(y)  are  found  by  equating  the  derivative  to  zero: 


dP(y)  = dP(y)  diL.  = o 

dy  du  dy 


(H-27) 


resulting  in  minima  p(y)  -»0  for  y — eo  (u  = 0)  and  y — *+ «o  (u  — + «o  ),  and 
a maximum  at  the  most  likely  value: 


p(y/>  = "tr  = 1-694 


(n-28) 


i.  e. , for 


u(yf  ) = y/T  , 


Y,  = 0.15  . 


(11-29) 


(H-30) 


l 
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The  normalized  log-envelope  probability  density  function  p(y),  and  the  cor- 
responding false  alarm  probability  P (y  ) are  sketched  in  Figure  II- 3. 

X1  1 

It  may  seem  strange  that  the  most  likely  value  of  y does  not 
correspond  to  that  ox  u (u^  = 1;  y(u^  ) = 0)*  This  is  due  to  the  fact  that  the 
transformation  y = log  u is  non-iinear,  which  shows  particularly  clearly 
in  Equation  (II- 24). 

4.  False  Alarm  Rate  Control 

The  above  shows  that,  by  dividing  the  envelope  by  a running 
estimate  of  the  RMS  noise  level,  we  obtain  a stationary  envelope  detection 
statistic  R/RMS.  For  any  desired  FAR  we  then  can  set  an  envelope  thres- 
hold R^,  at  the  corresponding  level  above  the  RMS  noise  value.  Table  II- 1 
lists  the  FAR  values  computed  for  the  righthand  tail  of  Figure  II- 3. 

The  table  shows  that  for  envelopes  of  a Gaussian  noise  process, 
relatively  low  FAR  values  are  obtained  for  threshold  settings  of  12  dB  or 
more  above  the  RMS  noise  level. 


Our  recently  designed  envelope  detector  (Unger,  1978)  requires 
that  two  detection  criteria  be  met  simultaneously: 

• The  signal-plus-noise  envelope  in  a test  gate  of  given  length 
must  exceed  the  maximum,  lagging  noise  envelope  a given 
number  of  times  (e.  g.  , 10%  to  30%  of  the  time  in  a 4- sec  test 
gate) 

• The  ratio  between  the  first  'signal'  envelope  peak  in  a test  gate 
and  the  maximum  noise  envelope  must  exceed  a given  threshold 
value  (e.  g.  , 2 or  3 dB). 

We  observe  from  Figure  II-3  that,  if  the  envelopes  are  Rayleigh-distributed, 
the  maximum  noise  envelope  is  most  likely  10  to  12  dB  above  the  RMS  noise 


II- 10 


y 


1 


US  ABOVE  RMS 


L06CR/RMSD 


jQi 

) assuming  minimum  interval  between  independent  samples  = 0.  36  sec 

FIGURE  II -3 

NORMALIZED  LOG-ENVELOPE  DISTRIBUTION  AND  PROBABILITY  OF 
FALSE  ALARM  (SINGLE  DETECTION  TRIAL)  FOR  GAUSSIAN 

NOISE  MODEL 
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TABLE  II- 1 

LOW  FAR  VALUES  FOR  SINGLE  DETECTION  TRIALS 


R/RMS 

(dB) 

PF 

FAR*) 

(FA/H) 

16 

- 8 

2.  00  x 10 

-5 

-4 

2x10 

14 

2. 00  x 10 

-3 

0.  2 

12 

1.34  x 10 

-2 

13.4 

10 

1. 69  x 10 

-2 

169 

8 

7.  88  x 10 

788 

6 

2.  04  x 10"1 

2040 

) Assuming  independent- sample  interval  = 0.36  sec. 


I 

] 

H-12  I 


value.  A single  detection  trial,  i.  e. , requiring  that  the  maximum  noise  en- 
velope is  exceeded  only  once,  then  would  give  a FAR  of  169  to  13.  4 FA/H. 
Multiple  detection  trials,  together  with  the  second  criterion,  will  reduce  the 
FAR.  In  an  initial  evaluation  of  this  envelope  detector  on  NORSAR  single- site 
data,  varying  the  threshold  settings  as  indicated  above,  i.  e.  , 10%  to  30%  for  the 
first  criterion,  and  2 to  3 dB  for  the  second  criterion,  yielded  FARs  of  7 to  20 
FA/H.  As  will  be  described  in  Section  III,  applying  the  envelope  detector  with 
thresholds  of  30%  for  the  first  criterion  and  3 dB  for  the  second  criterion  to 
one  hour  of  KSRS  single -site  noise  data  resulted  in  one  false  alarm.  Consider- 
ably more  noise  and  signal  data  need  to  be  processed  to  more  accurately  assess 
this  detector's  operating  characteristics.  This  is  discussed  further  in  Section 
IV. 


E.  FAR  CONTROL  WITH  THE  INSTANTANEOUS  NOISE  POWER 

1.  Normalized  Instantaneous-Noise  Distribution 

Below,  we  will  consider  an  alternative  detection  statistic,  the 
2 

instantaneous  power  n (t)  of  a zero-mean,  Gaussian  noise  process  n(t)  with  total 
2 

power  (r^  = N.  As  described  in  the  first  part  of  this  section,  the  r.  v.  n(t)  can 
be  normalized  by  substituting  v = n/ff^  to  obtain  a stationary  N(0,  1)  distribution: 


P(v) 


(II-  31) 


2.  Normalized  Log-Instantaneous-Power  Distribution 


We  now  will  study  the  distribution  of  the  base-ten  logarithm  of 

2 

the  normalized  instantaneous  power  v (t),  and  make  the  transformation 


w 


log  v 


(II-  32) 
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so  that 


dw  _ 2k 
dv  v 

where 

k = log  e = 0.  43429  as  defined  previously. 


(11-33) 


According  to  the  fundamental  theorem  (Papoulis,  1965;  p.  126),  we  find  the 
probability  density  function  for  w by 


n 


p(w)  = 2 

i=l 


P(vt) 

l8’(vi>! 


where 


g(v)  = w = log  v 

dw  2k 


g'(v)  = 


dv 


v are  the  roots  of  g(v),  i = 1 , 2,  ....  n. 


(11-34) 


In  our  case  the  roots  are 


so  that 


w 

2 


p(w) 


2k  y/In 


(II-35a) 
(II- 3 5b) 


+ 


2k  y/Tn 


e 


2 

v 

2 

2 

(11-36 ) 
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~'/T  ’ . 


I V1  I = I V2  I " V1  * 


-1  w 


(H- 37) 


p,”>  = ‘WSV  e 


(U-38) 


The  maxima  and  minima  of  p(w)  are  found  by  equating  its  derivative  to  zero: 


dP(w) 

dw 


= 0 , 


(11-39) 


resulting  in  minima  p(w)  —0  for  w — oo  (v  = 0)  and  for  w — -+«  (v  -►  + « ), 


and  a maximum 


p<w/ 1 = 


(11-40  ) 


at  the  most  likely  value 


wi  = 0 


(H-41  ) 


v(w#  ) = 1 


(U-42  ) 


The  functions  p(w)  and  P (w)  are  sketched  in  Figure  II-4.  Like  the  log- 

£ 

envelope  distribution,  the  log-instantaneous-power  distribution  is  non- symmet- 
ric (skewed),  and  its  most  likely  value  is  shifted  relative  to  that  of  v (vj  = 0; 
w(vj  ) — ► - co  ),  due  to  the  non-linear  logarithmic  transformation. 

3.  False  Alarm  Rate  Control 

Figure  II-4  shows  that  for  a Gaussian  noise  process  a desired 
FAR  is  obtained  by  controlling  the  instantaneous-power  threshold  relative  to 
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INST  PWR  is  the  instantaneous  power 
LTA  PWR  is  the  long-term  average  power 

*)  Assuming  minimum  interval  between  independent  samples  = 0.36  sec 

FIGURE  II -4 

NORMALIZED  LOG-INSTANTANEOUS-POWER  DISTRIBUTION  AND  PROBABILITY 
OF  FALSE  ALARM  (SINGLE  DETECTION  TRIAL)  FOR  GAUSSIAN  NOISE  MODEL 
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the  average  noise  power.  For  a given  threshold  wT»  the  FAR  is  predictable 
and  constant.  Low  FAR's  ( < 10  FA/H)  are  obtained  for  wT  = 1.  1,  corres- 
ponding to  instantaneous -power  thresholds  of  11  dB  or  more  above  the  average 
noise  power. 

Comparing  the  normalized  log-instantaneous-power  distribution 
with  the  normalized  log- envelope  distribution,  we  note  that  the  variance  of  the 
latter  is  much  smaller.  This  is  to  be  expected,  since  the  envelope  of  a wave- 
form, of  course,  has  less  variation  than  its  instantaneous  value.  This  char- 
acteristic could  make  the  envelope  a preferred  detection  statistic.  We  further 
observe  that  the  tail  of  the  normalized  log-instantaneous-power  distribution 
rises  less  sharply  with  decreasing  values  than  does  the  tail  of  the  log-envelope 
distribution.  At  first  thought,  one  might  interpret  this  as  though  the  instantane- 
ous power  would  provide  a more  favorable  FAR  control.  However,  we  do  not 
know  the  parameter  distributions  for  signal- plus -noise.  Determination  of  the 
difference  in  operating  characteristics  for  these  two  detection  statistics,  there- 
fore, requires  further  analysis  and  probably  empirical  evaluation. 

We  have  demonstrated  that,  for  Gaussian  noise,  control  of  the 
FAR  can  be  effectively  obtained  with  instantaneous  detection  parameters  derived 
from  a Gaussian  process.  It  is  not  necessary  that  a detection  parameter  it- 
self is  a Gaussian  r.  v.  Therefore,  for  Gaussian  noise,  averaging  of  an  instan- 
taneous detection  parameter  such  as  the  envelope  or  the  instantaneous  power, 
with  the  purpose  of  obtaining  a near -Gaussian  detection  statistic,  is  not  manda- 
tory in  the  design  of  a controlled  FAR  detector.  Averaging,  however,  may  en- 
hance the  performance  of  a detector  in  other  ways.  This  will  be  discussed  in 
Section  IV. 
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SECTION  in 
EMPIRICAL  DATA 


A.  INTRODUCTION 

In  this  section  we  discuss  detection  statistic  distributions 
compiled  from  empirical  data.  First,  we  show  log-envelope  distributions 
obtained  from  one  hour  of  KSRS  single-site  SP  noise  data,  and  discuss  the 
goodness -of-fit  relative  to  the  theoretical  distribution  derived  in  Section  U. 
These  data  are  then  augmented  with  observations  previously  reported.  Next, 
we  consider  power  data  obtained  by  Swindell  and  Snell  (1977)  in  the  course  of 
designing  and  evaluating  an  automatic  power  detector. 


B.  EMPIRICAL  ENVELOPE  DISTRIBUTIONS 
1.  Histogram  Generation 

The  previously  mentioned  envelope  detector  (Unger,  1978)  was 
adapted  to  handle  consecutive  records  as  uninterrupted  data.  In  this  form,  it 
was  applied  to  one  hour  of  uninterrupted,  unfiltered,  KSRS  single- site,  SP 
noise  data,  using  the  operating  parameters  listed  in  Table  IU-1.  These  pa- 
rameter values  are  similar  to  those  used  in  the  initial  envelope  detector  eva- 
luation, as  reported  in  the  above  mentioned  reference.  The  parameter  values 
listed  in  Table  III-l  are  also  in  use  at  the  present  time  in  the  evaluation  of  a 
multivariate  seismic  discrimination  package  (Sax  et  al.  , 1978).  The  one- 
hour  KSRS  noise  test  resulted  in  one  false  alarm. 

The  program  was  also  adapted  to  produce  histograms,  for 
102.  4-sec  data  segments,  of  the  log-envelope,  the  instantaneous  frequency 
and  the  mean  frequency,  using  various  sample -interval  lengths.  Below,  we 
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TABLE  ttl-l 

ENVELOPE  DETECTOR  OPERATING  PARAMETERS 


Warm-up  time  : 

60  seconds 

Peak-noise  tapering  time  constant  : 

180  seconds 

Peak-signal-to-peak-noise  threshold  : 

3 dB 

Signal  test  gate  : 

4 seconds 

Required  relative  frequency  of  detections  in  test  gate: 

30% 

Sample  interval  : 

0. 1 second 

will  analyze  the  log-envelope  histograms  obtained  from  sampling  consecutive 
data  segments  every  0.4  sec.  The  combination  of  KSRS  noise  spectrum  and  in- 
strument response  (Prahl  et  al. , 1975)  suggests  that  the  effective  bandwidth  of 
the  KSRS  noise  seismograms  approximately  equals  2.  5 Hz.  A sampling  interval 
of  0.4  sec  then  should  give  independent  samples  (Schwartz  et  al. , 1966).  This 
results  in  256  independent  envelope  samples  per  histogram.  The  transfer  of  de- 
tection parameters  from  one  record  to  the  next  requires  an  overlap  equal  to  the 
length  of  the  signal  test  gate,  in  this  case  4 sec.  The  start  times  of  consecutive 
histograms,  therefore,  are  98.4  sec  apart  rather  than  the  full  102.4  sec. 

2.  Histogram  Presentation 

The  resulting  histograms,  reflecting  a total  of  one  hour  of  noise 
data,  are  presented  in  Figure  III—  1 . The  log-envelope  values  have  been  norm- 
alized, as  described  in  Section  II,  by  subtracting  the  base-ten  logarithm  of  the 
RMS  noise  level  measured  over  the  entire  102.4-sec  data  segment.  The  upper 
axis  shows  the  corresponding  noise  envelope  values  in  dB  above  the  RMS  noise 
level.  The  histogram  at  11:52:12.0  (Figure  HI-1)  contains  the  only  • alarm'  re- 
ported by  the  detector.  The  corresponding  data  segment  is  shown  in  Figure 
IH- 2. 


Also  plotted  in  each  histogram,  as  a solid  curve,  are  the  rela- 
tive frequencies  found  from  the  theoretical,  normalized  log-envelope  distribu- 
tion for  a Gaussian  noise  process: 

Pj  = Ayp(y.)  , (HI-1) 

where 

i is  the  bin  index, 

Yj  is  the  value  at  the  center  of  bin  i , 

Ay  is  the  bin  width  (here  = 0.  1 ), 

p(Yi)  i*  the  probability  density  value  for  y^  as  given  by  Equa- 
tions (H-25)  and  (U-26),  and  in  Figure  II-3. 
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FIGURE  III- 1 

NORMALIZED  LOG-ENVELOPE  DISTRIBUTIONS  FOR  SP  NOISE,  KSRS  SITE  1,  9 JUNE  1976 
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Because  of  the  normalization,  this  theoretical  curve  is  identical  in  all  histo- 
grams. 

3.  Goodness -of- Fit  Testing 


Four  distribution  test  values  are  annotated  in  each  histogram. 

First,  the  value  SIG/MU  is  the  ratio  cr  / fi  , i.  e.  , the  ratio  of  envelope 

R R 

standard  deviation  over  envelope  mean,  computed  over  the  entire  102.  4- sec 
data  segment.  As  shown  in  Section  II,  for  a true  Rayleigh  distribution,  this 
ratio  should  equal  0.  525,  independent  of  the  noise  power. 


A second  test  is  the  chi-square  goodness -of- fit  test,  as  de- 
scribed in  standard  textbooks  on  statistics  (e.  g.  , Lindgren  et  al.  , 1978): 


£ <fi  - np/ 

i=l  1 1 


npi 


(HI- 2) 


where 

k is  the  number  of  bins  used  in  the  chi-square  test, 

f^  is  the  measured  frequency  of  occurrence  in  bin  i , 

n is  the  total  number  of  samples  used  to  compose  the  histo- 
gram, and 

p.  is  the  hypothetical  relative  frequency  as  given  by  Equa- 
tion (ni-1). 

2 

To  avoid  that  the  X -values  become  inflated  by  low  frequency-of-occurrence 
values  in  the  tails  of  the  distribution,  we  only  use  bins  in  the  interval  where 
the  theoretical  frequency  of  occurrence  is  greater  than  five,  as  suggested  in 
the  literature;  i.  e.  , for  the  values 

-0.5  4 y 4 0.4  , (IH-3) 

so  that  k=l 0 (9  degrees  of  freedom).  Disadvantages  of  the  chi-square  test  are: 
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• The  resulting  X distribution  depends  on  the  distribution  func- 
tion being  tested 

• Only  part  of  the  distribution  is  tested. 

Especially  because  of  the  second  item,  the  chi-square  test  is  not  the  best  one 

to  test  a continuous  distribution  such  as  our  theoretical  log-envelope  distribu- 
2 

tion.  Yet  the  X -values  were  easily  computed  and,  since  they  should  indicate 
at  least  a relative  goodness -of- fit,  they  were  taken  along,  and  annotated  as 
CHI  SQ  values  in  the  histograms. 

A better  test  for  continuous  distributions,  according  to  the 


literature,  is  the  Kolmogorov- Smirnov  goodness-of-fit  test: 


D s max  I F (y)  - F (y) 
n y I n o’ 


(HI-4) 


where 


F (y)  is  the  measured  cumulative  distribution  function  of  the 
n 

r.  v.  y , and 

F (y)  is  the  cumulative  distribution  function  to  be  tested,  i.  e.  , 
o’ 

the  hypothetical  cumulative  distribution  function. 

Thus,  equals  the  maximum  absolute  deviation  between  the  two  cumulative 
distribution  functions.  This  test  has  the  advantage,  that  it  automatically  cor- 
rects for  the  case  where  the  frequency-of-occurrence  may  not  be  smoothly 
distributed  over  the  chosen  bins,  i.  e.  , if  the  frequency  in  one  bin  is  too  high, 
but  correspondingly  low  in  the  next  bin,  the  Kolmogorov- Smirnov  test  still 
will  reflect  a 'good'  fit,  whereas  the  chi-square  test  will  not.  Another  im- 
portant advantage  is  the  fact  that  it  is  independent  of  the  distribution  being 
tested,  and  is  only  a function  of  the  sample  size  n . The  Kolmogorov- 
Smirnov  test  values  are  annotated  as  KOL>-SM  values  in  the  histograms. 

Finally,  we  recall  our  observation  in  Section  II,  that  the  ratio 
between  the  maximum  envelope  value  A and  the  RMS  noise  value  most  likely 
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is  10  to  12  dB  for  Gaussian  noise.  This  ratio  is  also  annotated  in  each  histogram. 


The  values  of  the  above  test  parameters  are  listed  in  Table  III- 2; 
their  distributions  are  plotted  in  Figure  III- 3.  The  apparent  relations  among 
these  various  test  parameters  are  shown  in  Figure  HI-4,  in  which  the  values  of 
SIG/MU,  CHI  SQ  and  A/RMS,  respectively,  are  plotted  against  the 
corresponding  KOL-SM  values.  The  ideal  SIG/MU  value  for  a Rayleigh  distri- 
bution of  envelopes,  and  the  CHI  SQ  and  KOL-SM  0.05  and  0.20  significance 
levels  a are  indicated  in  both  figures.  The  only  'alarm'  reported  is  indicated 
with  an  asterisk. 

4.  Discussion 

Our  discussion  will  be  based  on  how  well  our  theoretical  distri- 
bution can  be  used  in  establishing  adequate  control  of  the  FAR.  We  will  not  dis- 
cuss the  power  of  the  test  for  retaining  the  hypothesis  that  seismic  noise  is  Gaus- 
sian. The  goodness -of-fit  test  parameters  are  given  as  references  for  judging 

the  feasibility  of  using  the  presumed  distribution  as  a model  for  FAR  control. 

As  stated  earlier,  the  feasibility  of  FAR  control  hinges  upon  the 
possibility  of  obtaining  a stationary  detection  statistic  whose  distribution  can  be 
described,  to  a good  approximation,  in  closed  form.  The  approximation  must 
be  particularly  good  in  the  right  hand  tail  of  the  distribution. 

From  visual  inspection  of  Figure  III- 1 , we  believe  that  the  nor- 
malized Gaussian  noise  model,  with  the  ensuing  normalized  Rayleigh  distribu- 
tion for  normalized  noise  envelopes,  is  a reasonably  good  fit  to  the  observed 
occurrences,  at  least  for  the  one-hour  KSRS  single-site  noise  sample  tested. 

The  normalization  process  consists  of  dividing  the  envelope  values  by  the  RMS 
value  of  the  zero-mean  noise  seismogram.  Thus,  to  establish  adequate  FAR 
control  for  a given  noise  frequency  band,  we  merely  need  to  keep  a running  est- 
imate of  the  RMS  noise  value;  it  does  not  require  averaging  of  the  instantaneous 
amplitude  (envelope)  values.  Averaging  will  change  the  detector's  operating 
characteristics  (Section  IV). 
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GOODNESS-OF-FIT  TEST  PARAMETER  VALUES 


contains  detector-: 


26.  1 (1  value) 


In  80%  of  the  cases,  the  hypothetical  distribution  correctly  de- 
scribes at  least  the  right  hand  tail  of  the  empirical  envelope  distributions;  there, 
the  theoretical  curve  goes  virtually  through  the  bin  centers  at  the  correspond- 
ing relative  frequency  of  occurrence.  Exceptions  are  listed  in  Table  III-3,  with 
the  test  parameters  repeated  from  Table  III- 2.  It  includes  an  evaluation  of  the 
difference  in  FAR  these  tail  discrepancies  could  cause  if  we  would  set  a thres- 
hold at  1 1 dB  above  RMS.  The  remarks  in  Table  III- 3 describe  the  major  fea- 
ture of  the  particular  tail  discrepancy.  For  thresholds  of  12  dB  or  higher  there 
would  not  be  an  increase  in  the  FAR,  except  at  12:06:57.6  (Figure  III- 1 ) where 
a value  of  12.5  dB  above  RMS  was  recorded.  The  value  of  0.4  extra  false 
alarms  listed  in  Table  III- 3 for  this  case  is  based  on  102.4  sec  of  data.  As  we 
see  in  Figure  III- 1 , no  other  histogram  shows  envelope  values  above  12  dB 
above  RMS  in  the  entire  one-hour  data  segment,  suggesting  a FAR  of  only  0.  4 
FA/H  for  a single  detection  trial  at  a threshold  of  12  dB  above  RMS. 

Considering  Figures  III- 3 and  III-4,  we  note  first  that  the  tab- 
ulated significance  levels  correspond  fairly  well  with  the  cumulative  relative 
frequencies,  i.  e.  , 15%  of  the  Kolmogorov-Smirnov  test  values,  and  22.  5%  of 
the  chi-square  test  values  are  greater  than  the  respective  20%  significance 

levels.  Second,  it  seems  that  the  SIG/MU  values  start  to  diverge  beyond  the 

2 

KOL-SM  20%  significance  level.  Third,  besides  the  extreme  value  X = 26.  1, 
only  one  other  test  value  falls  beyond  the  KOL-SM/  CHI  SQ  20%  joint  signifi- 
cance level.  Fourth,  there  seems  to  be  adequate  correlation  between  CHI  SQ 
and  KOL-SM  test  values,  but  the  spread  in  this  relationship  increases  consid- 
erably with  increasing  test  values.  Finally,  we  observe  that  of  the  envelope 
distributions  with  tail  discrepancies  listed  in  Table  III- 3,  either  the  CHI  SQ  or 
the  KOL-SM  test  values  are  near  or  above  the  20%  significance  level. 

The  waveforms  and  the  envelope,  phase  and  frequency  time 
traces  of  two  data  segments,  one  with  a very  good  fit  and  another  with  a poor 
fit,  are  displayed  in  Figure  III- 5.  Besides  some  difference  in  waveform  character 
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TABLE  III -3 

ENVELOPE  DISTRIBUTIONS  WITH  DEVIATING  TAIL 


Time 

SIG/MU 

CHISQ 

* 

AFA  ) 

Remarks 

11:34:09.  6 

0.  588 

14.  3 

0.075 

11.7 

0.  9 

bin -quantization  noise 

11:45:38.  4 

0.  588 

9.9 

0.069 

10.9 

1.  8 

visually  non-Rayleigh 

11:50:33.  6 

0.614 

11.8 

0.085 

11.4 

2.  1 

visually  non-Rayleigh 

12:00:24.0 

0.  562 

10.  6 

0.036 

11.  3 

1.4 

bin -quantization  noise 

12:06:  57.6 

0.  599 

12.  3 

0.067 

12.  5 

large  A/RMS  value 

12:18:26.  4 

0.  565 

11.  3 

0.051 

11.  3 

B 

non -Rayleigh? 

12:29:55.  2 

0.496 

14.  1 

0.036 

8.9 

B 

FAR  overestimation 

Total  expected  FAR  deviation: 

s|e  a|e 

8.0  FA/H  = 20%  ) 

) Number  of  false  alarms  above  theoretical  number  (=  1.  1)  for  102.4- 
sec  data  segment,  assuming  0.  36- sec  independent- sample  interval  and 
envelope  threshold  of  11  dB  above  RMS  noise;  single  detection  trial. 

) Number  of  false  alarms,  for  one-hour  KSRS  noise  sample,  above  theo- 
retical FAR  (=  40  FA/H)  for  11  dB  threshold;  single  envelope  detection 
trials. 


INPUT 

12:31:33.  6 

'Good'  fit 
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FIGURE  ID- 5 

DATA  SEGMENTS  WITH  GOOD  AND  POOR  FIT,  RESPECTIVELY, 
TO  GAUSSIAN  NOISE  MODEL 
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there  is  little  to  suggest  that  the  poor-fit  waveform  contains  a signal.  The 
same  was  observed  for  all  data  segments  in  the  one-hour  KSRS  noise  sample. 
Neither  did  any  of  the  instantaneous  amplitude,  phase  and  frequency  traces 
display  any  anomalies;  all  looked  very  similar. 

In  summary,  we  conclude  that,  for  this  one-hour  KSRS  noise 
sample,  the  Gaussian  noise  model  is  an  adequate  model  to  use  as  a basis  for 
FAR  control  in  envelope  detection.  For  a desired  FAR,  we  only  need  to  set 
the  envelope  threshold  at  a corresponding  level  above  a running  estimate  of 
the  RMS  noise  level.  For  instance,  an  envelope  threshold  of  12  dB  above  the 
RMS  noise  should  yield  a FAR  of  13  FA/H  for  a single  detection  trial.  Employ- 
ing multiple  detection  trials  and  additional  detection  criteria  can  strongly  re- 
duce this  FAR.  We  will  discuss  this  further  in  Section  IV. 

These  analytical  results  should  be  tested  with  at  least  24  hours 
of  noise  data.  They  also  need  to  be  confirmed  with  noise  envelope  distribution 
analysis  on  noise  samples  recorded  at  different  times  of  day  and  year,  and  at 
different  stations.  Finally,  the  interval  for  independent  sampling  of  a waveform 
needs  to  be  determined  as  a function  of  the  noise  frequency  band  (Section  IV). 

C.  OTHER  OBSERVATIONS 

In  Section  II  we  noted  that  the  most  likely  ratio  of  maximum- 
envelope -over -RMS  for  Gaussian  noise  is  10  to  12  dB.  D.  G.  Lambert  ob- 
served that  in  noise  waveforms  recorded  at  stations  of  the  Very  Long  Period 
Experiment  (VLPE),  the  ratio  between  maximum  noise  amplitudes  and  RMS 
noise  values  was  rather  constant  (Lambert  et  al.  , 1973).  This  observation 
was  confirmed  for  several  VLPE  stations  when  studying  the  feasibility  to  de- 
termine a station's  detection  capability  from  noise  (Unger,  1974).  See  Table 
m-4.  Since  the  maximum  instantaneous  noise  value  in  general  equals  the 
maximum  noise  envelope  value  within  about  one  dB,  this  suggests  that  those 
LP  waveforms  may  also  have  Rayleigh-distributed  envelopes  and,  therefore. 
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TABLE  IU-4 

AVERAGE  MAXIMUM-AMPLITUDE-OVER-RMS  RATIOS  FOR  VLPE  NOISE*) 


Station 

Aero- 

nym 

AA/RMS 

s.  d. 

A/ RMS 
(dB) 

No.  Noise 
samples 

Charter  Towers,  Australia 

CTA 

8.  61 

2.  14 

12.  7 

26 

Chiang  Mai,  Thailand 

CHG 

7.  34 

1.20 

11.  3 

31 

Kongsberg,  Norway 

KON 

6.  55 

0.  58 

10.  3 

32 

Kipapa,  Hawaii 

KIP 

6.73 

0.  80 

10.  5 

48 

Albuquerque,  New  Mexico 

ALQ 

6.  26 

0.  36 

9.9 

30 

La  Paz,  Bolivia 

ZLP 

6.  60 

0.  68 

10.  4 

21 

) Unger  (1974),  p.  HI-ll 
K 

) AA  = largest  peak-to-peak  amplitude  in  one-hour  noise  sample 
A = largest  peak  amplitude  in  one-hour  noise  sample 
s.d,  = standard  deviation  of  AA/RMS  measurements 


would  stem  from  a truly  Gaussian  noise  process.  Although  all  A/RMS  values 
in  Table  III-4  agree  with  those  for  Gaussian  noise,  those  for  CTA  and  CHG 
are  larger  and  have  considerably  greater  variances  than  those  for  KON,  KIP, 
ALQ,  and  ZLP.  This  may  suggest  that  noise  at  the  latter  stations  is  more 
closely  Gaussian  than  that  at  CTA  and  CHG. 

We  furthermore  remark  that  in  many  cases,  based  on  empirical 
observations,  a 10  dB  peak-over-RMS  ratio  is  accepted  by  analysts  as  a stand- 
ard value  for  noise  waveforms. 

D.  POWER  DATA 

Since  our  main  interest  in  this  study  is  in  envelope  statistics, 
no  special  effort  was  made  to  collect  distributions  of  the  instantaneous  power. 
However,  Swindell  and  Snell  (1977)  collected  data  for  short-term  averages 
(STA)  of  the  instantaneous  power  in  order  to  obtain  a Gaussian  detection  sta- 
tistic. Some  of  their  results  are  shown  in  Figure  HI- 6,  presenting  distribu- 
tions for  logSTA  - logLTA  values  obtained  with  0.  8,  1.  6,  and  3.  2- sec  STA 
averaging  times.  LogLTA  is  the  long-term  (1  to  2 min)  mean  of  the  logSTA 

values,  so  that  LTA  approximately  equals  the  long-term  average  of  the  noise 

2 

power,  and  is  comparable  to  the  total  power  N or  RMS  in  our  previous 
analyses.  The  distributions  were  compiled  from  two  hours  of  uninterrupted 
KSRS  beam  noise;  in  our  present  study  we  used  the  first  hour  of  this  same 
period  to  compile  our  KSRS  single- site  envelope  distributions. 

As  observed  by  the  authors,  the  distribution  is  more  skewed 
for  the  shorter  integration  times.  The  distribution  for  the  0.8-sec  integration 
time  takes  on  more  the  shape  of  the  instantaneous  power  distribution  given  in 
Figure  II- 6.  Note  in  particular  the  fact  that  the  most  likely  value  occurs  near 
0 dB,  consistent  with  our  derivations.  This  helps  to  explain  why  a zero-mean, 
near-Gaussian  statistic  is  obtained  by  averaging  just  a few  independent  samples 
of  instantaneous  power  (four  0.4-sec  samples  for  the  1.  6-sec  gate;  eight  for  the 
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DISTRIBUTIONS  OF  LOG(STA)  - LOG(LTA),  FOR  VARIOUS 
AVERAGING  TIMES.  KSRS  90°  BEAM.  9 JUNE  1976,  11:30-13:30 
(SWINDELL  AND  SNELL,  1977;  p.  III-10) 
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( 

3.  2- sec  gate).  As  reported  by  the  authors,  the  means  of  these  distributions 
occur  slightly  different  from  0 dB  due  to  the  fact  that  the  instantaneous  power 
is  first  averaged  before  the  base-ten  logarithm  is  taken,  rather  than  averag- 
ing the  logarithms  of  the  instantaneous  power.  This  difference  approximately 

2 

equals  0.215  (q-  / u ) , where  <r  and  fi  are  the  standard  deviation  and 

P P P P 

the  mean  of  the  instantaneous  power  (Unger,  1974;  Appendix  C). 

The  consistency  of  the  logSTA  - logLTA  distributions  over 
various  noise  samples,  as  reported  by  Swindell  and  Snell  (1977),  strongly 
suggests  that  the  Gaussian  noise  model  may  in  general  be  an  adequate  basis 
for  the  design  of  a detector  with  a controllable  FAR. 
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SECTION  IV 

DETECTOR  DESIGN  CONSIDERATIONS 
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A.  INTRODUCTION 

In  the  previous  two  sections  we  have  shown  that  the  Gaussian 
noise  model,  which  generates  waveforms  whose  envelopes  follow  the  Rayleigh 
distribution,  provides  a stationary  detection  statistic  for  which  adequate  con- 
trol of  the  FAR  can  be  derived.  Below,  we  will  discuss  briefly  a number  of 
ways  to  best  use  the  information  obtained.  The  underlying  feature  of  all  de- 
sign strategies  is  the  assumption  that  a certain  tolerable  false  alarm  rate  will 
be  imposed  through  a surveillance  system's  threshold  control  (Sax  et  al. , 1974; 
Unger  et  al.  , 1974).  For  the  given  FAR,  the  probability  of  detection  then  needs 
to  be  maximized.  This  is  basically  the  Neyman- Pear  son  decision  rule,  treat- 
ed in  the  classical  decision  theory  and  applied  by  various  authors  to  the  field 
of  communications  technology  (Selin,  1965;  Schwartz  et  al.  , 1966;  Van  Trees, 
1968). 

To  discuss  these  design  strategies  in  great  detail  is  beyond  the 
scope  of  our  present  task.  Our  discussion  merely  serves  to  suggest  various 
approaches  that  may  be  taken,  and  to  indicate  areas  requiring  further  research, 
in  the  process  of  optimizing  the  design  for  a front-end  SP  signal  detector.  For 
more  details  the  reader  is  referred  to  the  above  mentioned  literature. 

B.  SIGNAL  DETECTION  BY  GOODNESS -OF -FIT  TESTING 

The  results  obtained  in  Section  III  in  testing  the  goodness-of- 
fit  of  the  envelope  Rayleigh  distribution  suggest  that  the  test  values  SIG/MU, 

CHI  SQ,  KOL-SM,  and  possibly  A/RMS,  may  be  used  to  either  retain  or  reject 
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) the  null  hypothesis  of  the  Gaussian  noise  model.  Rejection  of  this  null  hypo- 

thesis at  a certain  significance  level  a , corresponding  to  a desired  FAR, 
then  would  imply  the  presence  of  a signal. 

► Such  a procedure  would  require  continuous  updating  of  the  norm- 

alized envelope  distribution.  The  normalization  process  requires  a running 
estimate  of  the  RMS  noise  level.  For  this  strategy  it  is  extremely  important 
to  determine  the  minimum  interval  between  independent  samples,  and  the  mini- 
mum number  of  samples  (the  shortest  possible  gate)  from  which  a valid  noise 
envelope  distribution  may  be  estimated.  Furthermore,  the  effect  of  signals  on 
the  envelope  distribution  must  be  investigated;  the  distribution  of  signal-plus- 
noise  envelopes  will  determine  the  decision  'power'  of  the  test,  i.  e.  , the  pro- 
bability of  detection,  P^.  In  particular,  we  would  like  to  know  if  a relatively 
weak  signal  would  effect  the  total  noise  distribution  to  an  extent  that  the  null 
hypothesis  can  be  rejected  at  a specified  significance  level. 

As  a final  remark  on  this  strategy,  it  must  be  realized  that  any 
form  of  sufficiently  non- Gaussian  noise  would  be  recognized  as  'signail'.  Such 
false  alarms  may  be  due  to  large  noise  fluctuation  from  cultural  sources, 
atmospheric  disturbances,  or  distributed  ocean  sources.  Given  a detection, 
the  analyst  must  perform  additional  waveform  analysis  to  determine  if  the  de- 
tection satisfies  the  definition  of  a 'signal'  as  it  is  to  be  processed  in  a sur- 
veillance system,  i.  e. , in  detection  association,  event  location  and  classi- 
fication, etc.  Too  many  'detections'  of  non-Gaussian  noise  could  cause  the 
FAR  to  exceed  a level  required  for  optimum  execution  of  the  above  seismic 

, network  functions. 

C.  SIGNAL  DETECTION  BY  LIKELIHOOD  RATIO  TESTING 

In  many  signal  detection  cases  the  optimum  test  is  a likelihood 
ratio  test  (LRT).  Let  p^(y)  and  Pq(y)  be  the  probability  density  functions  of  a 
r.  v.  y under  the  condition  that  a signal  is  present  or  absent,  respectively. 
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The  LRT  then  consists  of  forming  the  likelihood  ratio 

p^y) 

L<y>  * • ov-» 

comparing  L(y)  with  a threshold  K,  and  deciding 

Dg  (only  noise  present)  if  L<  K, 

Dj  (signal  + noise  present)  if  L > K. 

Under  the  Neyman-Pearson  decision  rule  the  choice  of  K depends  on  the 
tolerable  FAR. 

For  m multiple  observations  or  detection  trials  ^ = 
y^,  y^» ...»  ym»  the  likelihood  ratio  is  the  ratio  of  the  joint  probability  density 
functions: 
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If  the  observations  are  independent  this  becomes 

77  P^/j) 

L(X)  = —LL* . (IV- 3) 

/7P0(V 

i=l 

In  many  cases  the  distributions  p^  and  p^  contain  exponential 
factors,  some  of  which  may  be  common  to  both  distributions.  Cancelling  the 
common  factors,  and  taking  natural  logarithms  on  both  sides  of  Equation  (IV- 3) 
usually  reduces  the  likelihood  ratio  to  a relatively  simple  detection  statistic 
which  maximizes  For  given  = a • For  instance,  in  the  case  of  env- 

elope detection  of  a signal  with  constant  amplitude  in  Gaussian  noise  by  multiple 
observation,  likelihood  ratio  testing,  the  optimum  necessary  and  sufficient 
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detection  statistic  can  be  shown  to  be  the  mean- square  envelope  for  low  SNR 
and  the  envelope  mean  for  high  SNR  (Schwartz  et  al. , 1966).  This  test  is  uni- 
formly most  powerful  in  the  sense  that  it  is  independent  of  either  the  SNR  or 
the  signal  amplitude;  any  other  test  for  this  case  would  be  less  powerful.  The 
power  of  the  test,  i.e.  , increases  with  both,  the  SNR  and  the  number  of 

independent  observations.  One  can  also  say  that  the  test  is  optimum  in  the 
sense  that,  for  given  P and  P , it  minimizes  the  SNR  required  for  signal 
detection.  Typically,  every  doubling  of  the  number  of  independent  observations 
(up  to  a certain  limit)  lowers  the  required  SNR  by  about  3 dB  (Schwartz  et  al.  , 
1966). 

In  the  case  of  seismic  event  signal  detection,  however,  there 
are  various  factors  which  may  make  it  difficult  to  perform  an  optimum  test  in 
the  above  fashion.  First,  the  number  of  independent  observations  is  a function 
of  signal  duration  and  of  signal  frequency  which,  in  turn,  depend  on  type,  magni- 
tude, and  orientation  of  the  seismic  event  and  on  the  propagation  characteristics. 
Second,  in  general  we  deal  with  signals  which  have  certain  rise  and  decay  times. 
Also,  especially  in  the  case  of  most  nuclear  explosions,  there  is  constructive 
and  destructive  interference  caused  by  early  arriving  multiple  signals.  The 
signal  amplitude,  therefore,  is  far  from  constant,  and  optimizing  the  test  for 
this  case  is  a difficult  task,  requiring  further  analysis. 

These  problems  are  recognized  in  radar  technology  as  'fading'  . 
Several  techniques  have  been  developed  to  cope  with  this  problem  (Schwartz 
et  al.  , 1966;  part  III).  Among  these  are  diverse  space,  time  and  frequency 
sampling  of  the  signal  waveform.  These  techniques  also  hold  promise  in  seis- 
mic event  signal  detection.  For  instance,  since  fading  will  frequently  be  dif- 
ferent at  different  sites  of  a seismic  array,  appropriate  multichannel  signal 
processing  may  improve  the  detection  of  fading  signals.  In  time  diversity, 
repetitions  of  the  signal  transmission  serve  a similar  function.  In  frequency 
diversity,  finally,  a bank  of  narrowband  filters  could  be  employed,  and  outputs 
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combined,  under  the  assumption  that  the  fading  would  vary  in  different  frequency 
bands. 


D.  AN  ANALYST-TYPE  DETECTION  MODEL 

We  will  now  turn  to  the  envelope  detector  designed  recently 
(Unger,  1978)  and  attempt  to  evaluate  its  features  based  on  the  information  from 
Sections  n and  III.  This  detector,  designed  to  closely  follow  the  detection  logic 
usually  applied  by  an  analyst,  operates  on  two  detection  criteria  (also  mention- 
ed in  Section  II)  which  must  be  satisfied  simultaneously: 

• The  signal-plus-noise  envelope,  sampled  within  a test  gate  of  giv- 
en length,  must  exceed  the  maximum  lagging  noise  envelope  for  a 
specified  percentage  of  the  total  number  of  samples  (e.  g. , for  30%, 
or  twelve  0.  l-sec  samples,  in  a 4-sec  test  gate). 

• The  ratio  between  the  first  'signal-plus -noise'  envelope  peak  in 
the  test  gate  and  the  long-term  maximum  noise  envelope  must 
exceed  a given  value  (e.  g.  , 2 or  3 dB). 

Furthermore,  the  maximum  noise  envelope  test  value  is  forced  to  decay  with 
time.  This  enables  the  detector  to  follow  long-term  (1  to  2 min)  noise  level 
variations. 

The  expected  number  of  times  that,  in  a T-sec  test  gate,  in  the 
absence  of  signal,  the  envelope  exceeds  the  lagging,  long-term  maximum  noise 
envelope  A is 

M = ~ jT~  • Pf(A/RM§)  , (IV-4) 

where  4T  is  the  interval  between  independent  samples  and  P is  the  proba- 
bility  of  a false  alarm  of  a single  detection  trial,  as  derived  in  Section  II  for 
normalized  envelopes  of  a Gaussian  noise  process.  In  Table  IV- 1 we  list  a 
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TABLE  IV- 1 

EXPECTED  NUMBER  OF  TIMES  THAT  NOISE  ENVELOPE 
IN  4-SEC  TEST  GATE  EXCEEDS  LONG-TERM  MAXIMUM 
NOISE  ENVELOPE  A,  FOR  GIVEN  A/ RMS*) 


1 


A/RMS 

(dB) 

** 

PF  > 

*** 

M ) 

12 

1. 34  x 10"3 

1. 34  x 10'2 

10 

1.  69  x 10'2 

1.  69  x 10"1 

6 

7.  88  x 10~2 

7.  88  x 10"1 

) Assuming  0.4-sec  independent- sample  interval 

sfeafc 

) From  Table  II- 1 
) From  Equation  (IV-4) 
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number  of  M -values  as  a function  of  A/RMS  values,  assuming  T = 4 sec  and 

4T  = 0.4  sec,  and  using  the  P -values  listed  in  Table  II- 1.  In  evaluating 

F 

the  envelope  detector  (Unger,  1978)  we  used  threshold  values  of  10%  and  30% 
for  the  first  detection  criterion.  For  T = 4 sec  and  AT  = 0.4  sec,  this  cor- 
responds to  M = l and  M = 3,  respectively.  In  view  of  Table  IV- 1 these  val- 
ues seem  quite  reasonable.  Only  in  the  unlikely  case  that  high  noise  envelopes 
are  clustered  within  a narrow  time  interval  would  a false  alarm  be  generated. 
In  contrast,  for  signals,  we  expect  that  relatively  high  envelopes  are  sustained 
over  at  least  a relatively  short  interval.  It  is  this  contrast  which  enhances  the 
power  of  our  hypothesis  testing,  or,  equivalently,  the  probability  of  detection. 

We  will  now  focus  on  the  second  criterion.  The  probability  that 
a noise  peak  in  the  test  gate  exceeds  the  lagging  long-term  noise  peak  A by 
B dB  is 

PB  * 1 * (*-VW)T/4T  • <IV-5' 


where 

y = log  A/RMS  + B/20  . (IV-6) 

/» t D 

For  instance,  according  to  Table  IV-1,  if  A/RMS  = 10  dB,  a SNR  threshold 
of  B = 2 dB  would  result  in  a probability  of  0.  0132  that,  in  the  absence  of  sig- 
nal, the  long-term  noise  peak  is  exceeded  by  2 dB  by  the  noise  peak  envelope 
in  the  test  gate. 

We  realize  from  the  above  considerations  that  the  events  of 
noise  satisfying  each  criterion  individually  are  not  independent.  The  joint 
probability  of  noise  satisfying  both  criteria  simultaneously,  therefore,  must 
be  less  than  each  of  the  individual  probabilities  but  greater  than  their  product. 
Thus,  denoting  the  probability  of  noise  satisfying  the  first  criterion  by  Pc,  we 
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have 


PB  PC<PBC<min(PB'PC)  * 


(IV- 7) 


This  shows  a range  over  which  the  second  criterion  further  reduces  the  false 
alarm  rate  for  this  detector. 

Clearly,  for  a given  pair  of  threshold  settings,  the  false  alarm 
rate  of  our  detector  depends  on  the  fluctuation  of  noise  peak-to-RMS  ratios. 
Although  this  ratio  was  observed  to  be  relatively  constant,  its  fluctuation  oc- 
curs over  a range  where  the  probability  of  false  alarm  increases  rather  sharp- 
ly with  decreasing  A/RMS  as  evidenced  by  Tables  II- 1 and  IV- 1.  This  may 
cause  some  instability  in  the  detector's  FAR.  The  extent  of  this  instability  is 
difficult  to  assess  analytically;  based  in  part  on  preliminary  empirical  evalu- 
ation, the  FAR  will  probably  fluctuate  between  0. 1 and  10  FA/H  for  threshold 
settings  of  30%  and  3 dB  for  the  respective  detection  criteria. 

This  FAR  stability  may  be  judged  unsatisfactory  for  control  of 
a seismic  surveillance  system.  For  Gaussian  noise,  the  FAR  can  be  stabiliz- 
ed by  referencing  our  envelope  detection  criteria  to  the  RMS  noise  level  rather 
than  to  the  maximum  noise  envelope. 

Such  an  envelope  detection  algorithm  resembles  that  of  an  STA/ 
LTA  type  detector.  The  essential  differences  between  the  envelope  and  the 
STA/LTA  type  detectors  lie  in  the  choice,  the  normalization  procedure,  and 
the  further  employment  of  the  detection  statistic.  These  differences  are  sum- 
marized in  Table  IV- 2,  for  the  automatic  Seismic  Research  Observatory  (SRO) 
detector,  the  automatic  power  detector  designed  by  Swindell  and  Snell  (1977), 
and  the  above  suggested  envelope  detection  scheme.  The  latter  is  attractive  in 
the  sense  that,  in  the  case  of  Gaussian  noise,  and  employing  single  detection 
trials,  stable  FAR  control  requires  neither  short-term  averaging  of  the 
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instantaneous  detection  statistic,  nor  dividing  by  the  standard  deviation.  As 
pointed  out  previously,  employing  an  instantaneous  detection  statistic  enables 
multiple  detection  trials  within  a presumed  signal  gate,  which  can  considerably 
improve  the  detector  operating  characteristics. 

This  improvement,  however,  depends  on  the  number  of  independ- 
ent samples  available  in  the  test  gate.  This  number  is  determined  by  the  noise 
spectrum.  Variation  in  the  noise  spectrum,  therefore,  may  cause  some  in- 
stability in  the  FAR.  For  instance,  the  interval  between  independent  samples 
is  0.4  sec  for  a bandwidth  of  2.  5 Hz,  0.  5 sec  for  2.  0 Hz.  Thus,  for  a 4-sec 
test  gate,  a change  from  2.  5 to  2.  0 Hz  in  the  noise  bandwidth  would  change  the 
number  of  independent  trials  from  10  to  8,  and  increase  the  FAR  accordingly. 

It  may  be  feasible  to  obtain  a stable  FAR  for  this  envelope  detector,  by  mak- 
ing the  threshold  dependent  on  the  noise  bandwidth.  The  above  also  suggests 
that  the  use  of  multiple  filters  with  partially  overlapping  bandwidths  may  bene- 
fit the  detector's  operating  characteristics.  These  issues  need  further  evalu- 
ation. 


Further  improvements  of  an  analyst-type  detection  model  by 
expanding  the  number  of  detection  observables  were  proposed  previously 
(Unger,  1978).  These  observables  could  include  the  envelope  slope,  the  instan- 
taneous phase,  and  the  instantaneous  frequency.  The  distribution  of  these  para- 
meters under  noise-only  conditions  need  to  be  determined  before  the  feasibility 
of  their  incorporation  can  be  assessed. 

Together  with  the  initial  evaluation  results  obtained  with  the 
envelope  detector  in  its  present  stage,  the  above  indicates  that  it  should  be 
feasible  to  design,  if  necessary  by  slight  modification  of  the  present  algorithm, 
an  envelope  detector  with  a controllable  FAR,  operating,  for  instance,  at  1 
FA/H  for  approximately  80%  signal  detection.  This  projection,  of  course, 
must  be  confirmed  with  extended  empirical  evaluation.  Further  experiment- 
ation with  the  threshold  controls  may  lead  to  further  improvement  of  the 
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operating  characteristics.  For  fading  signals,  such  as  those  from  seismic 
events,  we  believe  that  the  present  envelope  detector  with  the  proposed  modi- 
fications will  have  near-optimum  operating  characteristics. 

As  seen  in  Section  III,  for  the  purpose  of  FAR  control,  the 
Gaussian  noise  hypothesis  can  be  retained  at  a significance  level  of  5%  or 
possibly  less.  Most  likely,  this  will  cover  a great  number  of  seismic  station 
noise  conditions;  this  needs  to  be  verified  empirically.  For  conditions  under 
which  the  Gaussian  noise  hypothesis  has  to  be  rejected  at  this  relatively  low 
significance  level,  the  stability  and  the  predictability  of  the  FAR  for  a given 
threshold  setting  may  become  questionable.  Even  under  these  conditions,  how- 
ever, we  expect  that  the  analyst-type  envelope  detector  model,  operating  with 
multiple  detection  trials  and  multiple  detection  criteria  relative  to  the  long- 
term maximum  envelope,  will  provide  close  to  the  best  possible  performance 
characteristics . 

E.  STA/LTA-TYPE  DETECTOR  MODELS 

Above,  we  mentioned  briefly  both  the  apparent  similarity  and 
the  essential  differences  between  the  envelope  detector  and  two  kinds  of  STA/ 
LTA-type  detectors.  The  experience  gained  with  employment  of  the  latter  may 
be  useful  in  the  design  of  an  envelope  detector.  In  our  previous  analyses,  we 
already  discussed  some  of  this  experience.  Below,  we  will  elaborate  on  the 
FAR  stability  of  STA /LTA-type  detectors. 

In  Section  II  we  found  that,  for  Gaussian  noise,  the  normalized 
instantaneous  power  and  the  normalized  instantaneous  amplitude  are  stationary 
detection  statistics.  The  same  can  possibly  be  derived  for  a normalized  ver- 
sion of  the  rectified  instantaneous  value.  The  normalization  process  consists 
of  dividing  the  instantaneous  detection  parameter  by  a long-term  average:  the 
long-term  average  noise  power  for  power  detection,  the  long-term  RMS  noise 
for  envelope  detection.  We  would  not  expect  that  taking  a short-term  average 
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would  change  the  stationary  character  of  the  normalized  detection  statistic. 

For  Gaussian  noise,  therefore,  we  would  expect  that  straight  STA/LTA  de- 
tectors have  a constant  FAR  for  a given  threshold  setting.  However,  straight 
STA/LTA  detectors  have  been  reported  to  suffer  from  FAR  instability  (Lacoss, 
1972;  Steinert  et  al. , 1975;  Swindell  and  Snell,  1977).  This  instability  must  be 
caused  by  either  of  two  conditions,  or  a combination  of  these: 

• The  ST  A is  non- stationary 

• The  noise  is  non-Gaussian. 

In  the  case  of  Gaussian  noise,  a non- stationary  STA  statistic 
may  arise  when  the  noise  spectrum  is  non- stationary;  the  number  of  independ- 
ent samples  of  an  instantaneous  detection  parameter  within  a waveform  gate  of 
given  length  then  would  vary  with  time.  Since  the  variance  of  the  STA  depends 
on  the  number  of  independent  samples  in  the  test  gate,  noise  spectrum  vari- 
ation would  lead  to  fluctuation  in  the  STA  variance.  This  makes  the  STA  a non- 
stationary detection  statistic,  causing  FAR  instability.  For  this  case,  the 
FAR  of  STA/LTA  detectors  can  be  stabilized  by  normalizing  the  STA  statistic 
not  only  relative  to  the  noise  power  level  (the  LTA),  but  also  with  respect  to 
the  STA  variance.  This  method  was  suggested  by  Lacoss  (1972),  and  imple- 
mented by  Swindell  and  Snell  (1977)  in  their  design  of  an  automatic  power  de- 
tector with  a constant  FAR. 

Since  averaging  a r.  v.  will  tend  to  make  it  Gaussian,  the  above 
process  may  also  work  well  in  the  case  of  non-Gaussian  noise.  In  particular, 
the  logarithm  of  the  STA  of  the  instantaneous  power  becomes  near-Gaussian 
for  relatively  short  averaging  times  (Lacoss,  1972;  Swindell  and  Snell,  1977; 
see  also  Sections  II  and  III  of  this  report). 

Finally,  we  remark  that  averaging  of  an  instantaneous  detection 
statistic  in  certain  cases  is  equivalent  to  performing  a likelihood  ratio  test, 
which  may  be  an  optimum  test  for  those  cases  (Section  II). 


The  difference  in  performance  between  a detector  employing 
short-term  averaging  with  proper  STA  normalization,  and  a detector  operating 
with  multiple  detection  trials  and  multiple  detection  criteria,  should  be  further 
investigated. 

F.  DETECTOR- REPORTED  ALARM 

The  envelope  detector  in  its  present  design  stage  (Unger,  1978), 
with  the  adaptations  and  operating  parameters  as  described  in  Section  III,  was 
applied  to  the  one-hour  KSRS  single- site  noise  sample  analyzed  in  this  study. 
The  detector  reported  one  alarm  for  this  one-hour  period.  Below  we  will  ana- 
lyze this  alarm  in  terms  of  the  test  parameters  described  in  Section  III. 

The  reported  alarm  occurred  at  11:53:16.  3, in  the  data  segment 
starting  at  11:52:12.0.  The  log-envelope  distribution  for  this  segment  is  shown 
on  page  4 of  Figure  III- 1 ; the  waveforms  are  presented  in  Figure  III-2.  The 
goodness -of- fit  test  parameters  for  this  segment  are  marked  with  an  asterisk 
in  Table  III- 2,  and  in  Figures  III- 3 and  III-4. 

We  note  that,  although  the  right  hand  tail  of  this  data  segment's 
envelope  distribution  fits  the  Rayleigh  distribution  very  well,  the  relative  fre- 
quencies to  the  left  of,  and  including  the  most  likely  value,  show  considerable 
deviation.  This  is  also  evidenced  by  the  extreme  SIG/MU  value,  the  high  KOL- 
SM  value,  and,  to  a lesser  extent,  by  the  CHI  SQ  value.  As  observed  for  the 
waveforms  of  other  data  segments  with  unfavorable  test  values,  however,  there 
is  little  indication  of  signal  presence  in  Figure  III-2,  suggesting  that  this  was 
a false  alarm. 
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SECTION  V 
SUMMARY 


Distributions  of  the  base-ten  logarithm  of  instantaneous  amp- 
litude (envelope)  values  from  one  hour  of  uninterrupted,  unfiltered  KSRS  sin- 
gle-site, SP  noise  were  analyzed.  The  one-hour  noise  sample  was  divided  into 
102.4-sec  data  segments.  Taking  0.4-sec  envelope  samples,  the  log-envelope 
distribution  was  compiled  for  each  data  segment,  in  2 dB  bins.  The  evaluation 
led  to  the  following  conclusions: 

• For  the  envelopes  of  this  one-hour  noise  sample  the  Rayleigh 
distribution,  which  is  the  theoretical  envelope  distribution  for 
waveforms  generated  by  a Gaussian  noise  process,  can  be  re- 
tained at  the  Kolmogorov-Smirnov  5%  significance  level. 

• Given  a Rayleigh  distribution,  a stationary,  normalized  detect- 
ion statistic  is  obtained  by  dividing  the  envelope  values  by  a 
long-term  (1  to  2 minutes)  running  estimate  of  the  RMS  noise 
distribution. 

• At  the  5%  significance  level  of  retention,  the  Rayleigh  distribu- 
tion seems  an  adequate  model  for  deriving  stable  control  of  the 
FAR  of  an  envelope  detector.  For  a single  detection  trial,  a 
tolerable  false  alarm  rate  can  be  maintained  merely  by  setting 
the  envelope  signal  detection  threshold  at  a corresponding  level 
above  the  RMS  noise  level.  For  instance,  a threshold  setting 

of  12  dB  above  RMS  should  yield  13  false  alarms  per  hour.  This 
is  yet  to  be  confirmed  empirically. 
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Employing  multiple,  independent  detection  trials  and  multiple  f 

detection  criteria  can  further  reduce  the  false  alarm  rate.  The 
amount  of  false  alarm  rate  reduction  depends  on  the  noise  fre- 
quency hand.  Therefore,  noise  frequency  variation  may  cause 
some  false  alarm  rate  instability.  The  optimum  operating  para- 
meters will  have  to  be  determined  empirically. 

For  Rayleigh-distributed  envelopes  in  an  approximately  1.  5-minute 
gate,  the  theoretically  most  likely  value  of  the  ratio  between 
the  maximum  noise  envelope  or  the  maximum  instantaneous  noise 
value,  and  the  RMS  noise  level  is  10  to  12  dB.  This  value  is 
routinely  observed  by  seismogram  analysts.  Since  for  Gaussian 
noise  this  ratio  is  nearly  constant,  a detector  operating  rela- 
tive to  the  long-term  noise  envelope  maximum,  e.  g.  , with  a 
threshold  of  3 dB  above  this  maximum,  may  have  a near-con- 
stant false  alarm  rate.  It  may  exhibit  a relatively  stable  per- 
formance also  under  non-Gaussian  noise  conditions. 


A recently  developed  envelope  detector  (Unger,  1978)  operating 
relative  to  the  long-term  noise  envelope  maximum,  and  using 
multiple  detection  trials  and  multiple  detection  criteria  closely 
following  an  analyst's  detection  logic,  reported  one  false  alarm 
over  the  one-hour  noise  sample.  Determination  of  reliable  op- 
erating characteristics  for  this  detector  model  requires  further 
empirical  evaluation. 

A detection  statistic  need  not  be  Gaussian  in  order  to  enable 
efficient  false  alarm  rate  control;  it  merely  has  to  be  stationary. 
Parameters  derived  from  a Gaussian  process  can  conveniently 
be  transformed  into  stationary  detection  statistics.  Examples 
are  normalization  of  the  Rayleigh-distributed  envelope  and 
normalization  of  the  instantaneous  power. 
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• In  view  of  the  above,  in  the  case  of  Gaussian  noise,  taking  a 
short-term  average  (STA)  of  an  instantaneous  detection  statis- 
tic is  not  mandatory  in  the  design  of  a controlled  false  alarm 
rate  detector.  In  fact,  fluctuations  in  the  noise  frequency  band 
would  cause  the  STA  statistic  to  be  non- stationary  (Lacoss,  1972; 
Steinert  et  al.  , 1975);  this  explains  the  reported  false  alarm 
rate  instability  in  straight  STA/LTA  detectors.  The  STA/LTA 
false  alarm  rate  can  be  stabilized  by  proper  normalization  of 

the  STA  statistic  (Lacoss,  1972;  Swindell  and  Snell,  1977). 

• Short-term  averaging,  followed  by  STA  normalization,  may  en- 
hance the  detector  operating  characteristics  in  the  sense  of  per- 
forming a likelihood  ratio  test  for  multiple  trials  or  observations 
in  the  detection  of  signals  of  constant  amplitude.  It  also  may  pro- 
vide a means  of  controlling  the  false  alarm  rate  for  non-Gaussian 
noise. 

• For  signals  of  non-constant  amplitude,  such  as  most  seismic  sig- 
nals, it  is  difficult  to  analytically  optimize  the  detector  design. 

A combination  of  analytical  and  empirical  evaluation  seems  neces- 
sary. In  particular,  the  strategy  of  normalized  STA/LTA  detec- 
tion (e.  g. , Swindell  and  Snell,  1977)  should  be  compared  to  that 
of  employing  multiple,  independent  detection  trials  and  multiple 
detection  criteria  based  on  instantaneous  detection  statistics 
(Unger,  1978). 

• For  all  possible  detector  design  strategies,  it  is  extremely  im- 
portant to  determine  the  shortest  interval  between  independent 
samples  as  a function  of  the  noise  frequency  band. 

• The  suggested  noise  model  should  be  tested  on  noise  data  re- 
corded at  different  times  of  day  and  year,  and  at  a variety  of 
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stations.  Sufficiently  long  detector  runs  should  be  carried  out 
to  establish  the  validity  of  the  noise  model  in  predicting  the  oc- 
currences of  noise  alarms  in  the  tail  of  the  distribution,  i.  e.  , 
at  false  alarm  rates  sufficiently  low  for  optimum  network  oper- 
ation. 
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